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Abstract 

Real filaments are not perfectly homogeneous. Most of them have 
various materials composition and shapes making their stiffnesses not 
constant along the arclength. We investigate the existence of circular 
and helical equilibrium solutions of an intrinsically straight rod with 
varying bending and twisting stiffnesses, within the framework of the 
Kirchhoff model. The planar ring equilibrium solution only exists for 
a rod with a given form of variation of the bending stiffness. We show 
that the well known circular helix is not an equilibrium solution of 
the static Kirchhoff equations for a rod with non constant bending 
stiffness. Our results may provide an explanation for the variation 
of the curvature seen in small closed DNAs immersed in a solution 
containing Zn^+, and in the DNA wrapped around a nucleosome. 

Keywords: Kirchhoff rod model, inhomogeneous rod, planar ring, generalized 
helices, planar elasticae, closed DNA 



1 Introduction 

Static and dynamics of one-dimensional structures have been the subject 
of intense research in Engineering, Physics and Biology. Such filamentary 



objects are present in both nature and man-made devices. The stabihty of 
sub-oceanic cables [I1I2|, the tridimensional structure of biomolecules and 
bacterial fibers ^ , and the mechanical properties of nanosprings El 13 IE] 
are few of many interesting problems investigated within elastic rod models. 

In most cases, the rod is considered as having constant or uniform stiff- 
ness, but non-uniform rods have been considered in the literature. It has been 
shown that nonhomogeneous Kirchhoff rods may present spatial chaos |^[TU]. 
Domokos and Holmes studied buckled states of non-uniform elasticae [TT] . 
Deviations from the helical structure of rods due to periodic variation of 
the Young's modulus were verified numerically by da Fonseca, Malta and 
de Aguiar ^2]- Homogeneous and nonhomogeneous rods subject to given 
boundary conditions were studied by da Fonseca and de Aguiar in J^l- The 
effects of a nonhomogeneous mass distribution in the dynamics of unstable 
closed rods have been analyzed by Fonseca and de Aguiar ^1]. Goriely and 
McMillen ^S] studied the dynamics of cracking whips ^^1; Kashimoto and 
Shiraishi studied twisting waves in inhomogeneous rods, and Neuringer 
and Elishakoff studied natural frequencies of inhomogeneous rods. 

It is well known that straight, circular and helical rods are equilibrium 
solutions of the Kirchhoff model [ini EHl EI] • These solutions have been used 
to study the stability and buckling phenomena of loop elasticae ^11221; the 
twist to writhe conversion phenomenon in twisted straight rods [23 , and the 
stability of DNA configurations with or without self-contact j21], to name 
some examples. 
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Our aim here is to investigate the existence of circular and hehcal equihb- 
rium solutions for an intrinsically straight, inextensible, unshearable, and in- 
homogeneous rod with circular cross- sect ion, within the Kirchhoff rod model. 

Establishing the stability (even in the case of simple equilibrium config- 
urations) of a rod is a non-trivial task (see, for example, Refs. [221 123 1211 
1221 12ni I2Z1 12H1)- Since the study of the stability analysis departs from the 
knowledge of equilibrium solutions, our work is the starting point for the 
study of the stability of equilibrium solutions of inhomogeneous rods. 

The inhomogeneity of the rod is considered through the bending and 
twisting coefficients varying along its arclength s, B{s) and C(s), respectively. 
We shall derive a set of non-linear differential equations for the curvature, 
kpi^s), and the torsion, tf{s), of the centerline of an inhomogeneous rod, and 
then impose the necessary conditions to finding circular/helical solutions. 
Circular (ring) solutions have kp = Constant and Tp = 0, while helical 
solutions have the curvature and torsion of the rod centerline satisfying the 
Lancret's theorem: kF{s)/TF{s) is constant [2nilSII|- The well known circular 
helix is the particular case where both kp and tf are constant. 

The fundamental theorem for space curves [20] states that the curvature 
and the torsion completely determine a space curve, but for its position 
in space. This is the reason for our choice of working with the curvature, 
kpis), and the torsion, tf{s), of the rod centerline. We shall show that the 
solutions for kp^s) and tf{s) depend only on the bending coefficient, B{s), 
an expected result since the centerline of the rod does not depend on the 
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twisting coefficient (see for instance, Neukirch and Henderson jSI])- 

In Sec. II we present the static Kirchhoff equations for an intrinsically 
straight rod with circular cross section and varying stiffness, and derive the 
non-linear differential equations for the curvature and torsion of the rod. In 
Sec. Ill we analyse the cases of null torsion (straight and planar rods). We 
show that the ring solution of the static Kirchhoff equations only exists for 
a particular form of variation of the bending stiffness. In the case of twisted 
rods with varying bending stiffness, this ring solution is the only possible 
planar solution of the Kirchhoff equations. In the case of non-twisted rods, 
there exist other types of equilibrium planar solutions that depend on the 
bending stiffness. In Sec. IV we use the Lancret's theorem for obtaining 
generalized helical solutions of the static Kirchhoff equations. We show that 
the circular helix is not an equilibrium solution of a rod that has non constant 
bending stiffness. We obtain generalized helical solutions that depend on 
the form of variation of the bending stiffness. As illustration, we compare 
the helical solution of a homogeneous rod with two types of helical solutions 
related to simple cases of inhomogeneous rods: (i) bending coefficient varying 
linearly, and (ii) bending varying periodically along the rod. In Sec. V we 
discuss the main results and possible applications to DNA. 
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2 The static Kirchhoff equations 

The statics and dynamics of long, thin, inextensible and unshearable elastic 
rods are described by the Kirchhoff rod model. In this model, the rod is 
divided in segments of infinitesimal thickness to which the Newton's second 
law for the linear and angular momenta are applied. We derive a set of 
partial differential equations (PDE) for the averaged forces and torques on 
each cross section and for a triad of vectors describing the shape of the rod. 
The set of PDE are completed with a linear constitutive relation between 
torque and twist. 

The central axis of the rod, hereafter called centerline, is represented by 
a space curve x parametrized by the arclength s. We commonly describe a 
physical filament using a local basis, {di, d2, ds}, which permits taking into 
account the twist deformation of the filament. This local basis is defined such 
that ds is the vector tangent to the centerhne of the rod (da = t = dx/ds), 
and di and d2 he on the cross section plane. The local basis is related to the 
Frenet frame {n, b, t} of the centerline through 



(di d2 d3) = (n b t) 



^ cos^ — sin^ ^ 

sin ^ cos ^ 
y ly 



(1) 



where the angle ^ is the amount of twisting of the local basis with respect to 
t. 
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We are interested in the equilibrium solutions of the Kirchhoff model, so 
our study departs from the static Kirchhoff equations For intrinsically 
straight isotropic rods, these equations are: 



F' = , (2) 
M' = F X dg , (3) 

M = B{s) h di + B{s) h d2 + C{s) h da , (4) 



where the following scaled variables were introduced: 



EI EI 1 

sL , F^F— , M->M— and h ^ h- ,i = 1,2,3 . (5) 

L L 



the vectors F and M are the resultant force and the corresponding moment 
with respect to the centerline of the rod, respectively, at a given cross section. 
The prime ' denotes differentiation with respect to s. ki are the components 
of the twist vector, k, that controls the variations of the director basis along 
the rod through the relation 

d^ = k X di , i = 1,2,3 . (6) 

ki and ^2 are related to the curvature of the centerline of the rod {kp = 



\/kl + k^) and k^ is the twist density. B[s) and C{s) are the bending and 
twisting coefficients of the rod, respectively. Writing the force F in the 
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director basis, 

F = /idi + /2d2 + /ads , (7) 

the equations (j2FEl) give the following differential equations for the compo- 
nents of the force and twist vector: 

f[-f2h + fsk2 = 0, (8) 

f2 + fih-f3ki = 0, (9) 

f3-flk2 + f2k, = 0, (10) 

(Bis) ki)' + {C{s) - Bis)) k2k,-f2 = 0, (11) 
{B{s) k^y - (Cis) - Bis)) k,ks + fi = 0, (12) 

(Cis) ks)' = . (13) 

The equation pHj) shows that the component M3 = C{s) k^ of the moment in 
the director basis (also called torsional moment), is constant along the rod, 
consequently the twist density k^ is inversely proportional to the twisting 
coefficient C{s) 

*s(.) = ^ . (14) 

In order to look for circular (ring) and helical solutions of the Eqs. (|5HT^ 
the components of the twist vector k are expressed as follows: 

ki = kpis) sin^ , (15) 
^2 = kp^s) cos^ , (16) 



k3 = + 



(17) 



where kpi^s) and Tp{s) are the curvature and torsion, respectively, of the rod 
centerhne, and ^ is given by Eq. ((T)). 

Substituting Eqs. (fT3HT7|) in Eqs. (jHHTSj). extracting fi and /2 from 
Eqs. flT^ and (fTTj) . respectively, differentiating them with respect to s, and 
substituting in Eqs. (jH)), © and (fTUj) . gives the following set of nonlinear 
differential equations: 

[M^kp - BkpTF]' - {EkpyTF = , (18) 
(B kp)" + A;^ r^[M3 - 5 r^] - /g A;^ = , (19) 
(EM'^F + /3 = 0, (20) 

where we have omitted the dependence on s to simplify the notation. Ap- 
pendix A presents the details of the derivation of Eqs. ()18fEn|l . 

The Eqs. (fTH^OUj) for the curvature, kp, and torsion, r^, do not depend 
on the twisting coefficient, C{s). Therefore, the centerline of an inhomo- 
geneous rod does not depend on the twisting coefficient unlike the case of 
homogeneous rods (see, for example, Eqs. (13) and (14) of Ref. [SI])- 

Langer and Singer [HBl have obtained a set of first-order ordinary differen- 
tial equations for the curvature and torsion of the centerline of a homogeneous 
rod that contains terms proportional to k'^p and t^. The Eqs. fjl8hOm) have 
the advantage of involving only terms linear m. kp and Tp. 
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3 Planar solutions of inhomogeneous rods 

Planar solutions are obtained by setting = in Eqs. (jl8H2(Jj) . The result 
is: 



M^k'p = 0, (21) 
{BkF)"-f,kF = 0, (22) 
(5M'^F + /3 = 0, (23) 

Eq. (PT|) implies two different sets of solutions: M3 7^ (twisted rods) 
and M3 = (non- twisted rods). Each case will be analysed separately. 



3.1 Twisted rod: M3 7^ 

As (M3 7^ 0), Eq. ()21|) leads to kpis) = Constant = K, which represents the 
circular solution (also known as planar ring solution). Eqs. (j22I) and ()23j] 
must be also satisfied, and their combination gives the following differential 
equation for the bending stiffness B{s): 

A. [B" + K^B] = , (24) 

yJjO 

or 

B" + K'^B = Constant = Cq , (25) 
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where we have set kp = K, and Cq is a constant of integration. Therefore, 
the twisted (M3 7^ 0) ring solution of an inhomogeneous rod only exists if 
its bending stiffness satisfies the differential equation (j25|) (forced harmonic 
oscillator). The bending stiffness (solution of ()25|) ) can be written in the 
form: 

B{s) = Ao cos{K s) + Bo sm{K s) + Cq/K'^ , (26) 

with Aq and Bq arbitrary constants. 

Therefore, only an intrinsically straight inhomogeneous twisted rod (M3 7^ 
0) with the bending stiffness given by Eq. (j^ . can display a planar ring 
configuration. 

3.2 Non- twisted rod: M3 = 

In this case, Eq. ()2H) is automatically satisfied so the curvature kpis) does 
not have to be constant. 

From Eq. (jl2I), we can obtain the following relation for fsis): 



Ms) - . (27) 



Differentiating with respect to s and substituting in Eq. (|^. we obtain 
the following non-linear differential equation for the curvature kpis): 

kFiBkp)'" - k'p{BkF)" + kliSkp)' = . (28) 
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If kp = Constant = K (the planar ring solution) the above equation gives 
the following differential equation for B: 



4- \B" + K^B] = , (29) 
as 



which is the same Eq. ()24|) and, therefore, gives the same condition for B{s) 
shown in Eq. In this case, fsis) is given by: 

/3(s) = -K^{AocosKs + Bo sin Ks) . (30) 

So, the non-twisted planar ring {kp = Constant) is also an equilibrium 
solution of the static Kirchhoff equations if the bending stiffness of the inho- 
mogeneous rod is given by Eq. ^E^ . 

Another solution of Eq. (j^Hj) for the curvature kpls) comes from making 
[Bkp)' = 0, which gives 



, , , Constant , , 



and, from Eq. (jZH), fsis) = 0. 

It is a very interesting solution because it relates the curvature of the rod 
centerline to the bending stiffness. Also, the component of the force tangent 
to the rod centerline, /a, is null for the set of solutions satisfying Eq. (p?T|l . 
In the particular case where B{s) is given by Eq. (pUj) . we can have two 
possible planar solutions: the planar ring /ci?(s) = K, and the solution for 
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which kpi^s) is given by Eq. (jHH). 

A kink is characterized by a piece of a filament having high curvature 
kp (or small radius of curvature, R = 1/kp). The Kirchhoff approach is 
not appropriate for modelling kinks in rods because in the derivation of the 
Kirchhoff equations it is assumed that the rod centerline radius of curva- 
ture is much larger than the cross-section radius pUj. However, Haijun and 
Zhong-can |33] used the Kirchhoff rod model to show that variations in the 
curvature of a closed rod may lead to instabilities responsible for the kink 
formation. Within the limits of validity of the Kirchhoff model, Eq. (13 1|) 
shows that pieces of a rod with small bending coefficient may lead to equi- 
librium configurations with large curvature. Therefore, our results can help 
explaining the onset of formation of kinks, since variations in the curvature 
of a filament can be originated by variations in the bending stiffness along it 
(Eq. (jHlll ). In Section V we comment on an application of this result to small 
closed DNAs. Figure ^ illustrates a planar equilibrium solution given by Eq. 
(p?T|l in the case of a rod with periodic variation of its bending stiffness. 

The straight rod {ki = 0, ^2 = 0) is a trivial solution of the static 
Kirchhoff equations, and is a particular case of planar solutions. In the 
inhomogeneous case, the twist density, fcs, of the straight rod is not constant. 
According to Eq. ()14|) . for a given external torque, M3, the twist is large for 
pieces of the rod where the twisting coefficient, C, is small. This result could 
help in determining the regions where plastic deformations or ruptures start 
to occur in rods subjected to large stresses. 
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Figure 1: A non twisted planar solution of the Kirchhoff equations for a rod 
with periodic variation of the bending stiffness. The arrows point to the 
regions where the curvature is increased. 

4 Helical solutions of inhomogeneous rods 

In this section, we are considering the equilibrium solutions for the rod cen- 
terline in which kp ^ ^ and rp 7^ 0. In order to find helical solutions for the 
static Kirchhoff equations, we apply the Lancret's theorem to the equations 
fll8fEn|) . We first write the Lancret's theorem in the form: 



where /? 7^ is a constant. Substituting Eq. (jHSI) in Eq. (fTH|) we obtain 



(32) 



t'f{M^ - B rp) - 2tf {B Tf)' = Q . 



(33) 
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Substituting Eq. (jH^ in Eq. (fT^ and extracting f^, we obtain 



f3=^-^^ + Tp{Ms-BTF) . (34) 



Differentiating with respect to s and substituting in Eq. (|2(Jj) we obtain 
the following differential equation for tf{s): 

(£-r „.^l,.,(B,,)-.0, (35) 

where Eq. was used to simplify the above equation. One immediate 

solution for this differential equation is 

{Btf)' = Q, (36) 

that substituted in Eq. (IHHj) gives 

t'p (Ms - Btf) = . (37) 

For non-constant r^?, the Eq. (|?fj) gives the following solution for Tp: 



rF{s) = —- . (38) 



If M3 = 0, the only solution satisfying the Lancret theorem is Tp{s) = 
which is the planar solution given in the previous section. Substituting Eqs. 
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and (jSHl) in Eq. ^ we obtain that 



fsis) = . 



(39) 



Substituting Eq. in (j22I) we obtain: 



(40) 



Substituting Eq. ((H) in Eq. 0, the force F becomes 



F = (/i cos^ - /s sin^) n + (/i sin^ + /a cos^ b + /a t 



(41) 



where {n, b,t} is the Frenet basis. Using the Eqs. ()57|) and (j3Hj) for fi and 
/2 (Appendix A), we obtain 



where in the inhomogeneous case /s must satisfy the Eq. ()2()|1 . 

Substituting Eqs. in the Eq. (021), it follows F = 0. Therefore, 

the helical solutions satisfying ()36|) are free standing. 

Now, we show that a circular helix cannot be a solution of the static 
Kirchhoff equations for an intrinsically straight rod with non-constant bend- 
ing stiffness. A circular helix has kp = Constant 7^ 0, and rp = Constant 



F 



(42) 
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7^ 0, in which case Eq. (fT5|) gives: 



2kFTFB' = . (43) 

Since kp ^ and Tp ^ 0, Eq. (j43|) will be satisfied only ii B' = 0, implying 
constant bending stiffness. Therefore, it is not possible to have a circular 
helix as a solution for a rod having non-constant bending coefficient. 

The solution for the curvature hp, Eq. (jlOI), and the torsion Tp, Eq. (jSH)), 
can be used to obtain the unit vectors of the Frenet frame {n, b, t} through 
the Frenet-Serret equations: 

t' = kpn , (44) 
n' = -kpt + Tph , (45) 
b' = -TfU . (46) 

By choosing the 2;-direction of the fixed cartesian basis as the direction of 
the helical axis, we can integrate t in order to obtain the three-dimensional 
configuration of the centerline of the rod. 

Figure 121 displays the helical solution of the static Kirchhoff equations for 
rods with bending coefficients given by 

FigEt: Bais) = 1, (47) 
FigEb: Bbis) = 1 + 0.007s , (48) 
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Fig|2t: B,{s) 



= 1 + 0.1 sm(0.04s + 2) 



(49) 



The case of constant bending (Eq. (|47|) ) produces the well known circular 
helix displayed in Fig. while Figs. |2bH2t display other helical structures 
obtained for rods with non-constant bending coefficients (Eqs. ()48H49|) ). 
These structures are clearly not circular helices but as they satisfy the Lancret 
theorem, Eq. ()32j) . they are called generalized helices. 



(a) (c) 




Figure 2: Helical solutions of the Kirchhoff equations using the Lancret's 
Theorem, (a) circular helix solution valid for a rod with constant bending 
coefficient Ba = 1 (@7j); (b) and (c) are generalized helices for inhomogeneous 
rod with bending coefficient given by Eqs. ()48|) and (jlHl), respectively. The 
parameters, in scaled units, are M3 = 0.05, and the total length of the rod 
is L = 130. kpiO) = 0.24 for the helical solutions displayed in panels (a) and 
(b), and ^^(0) = 0.22 for panel (c). 

The tridimensional helical configurations displayed in Fig. |21 were ob- 
tained by integrating the Frenet-Serret equations ()44ll4(jj) using the follow- 
ing initial conditions for the Frenet frame: t(s = 0) = (0, sin a, cos a), 
n(s = 0) = (—1,0,0) and b(s = 0) = (0, — cos a, sin a), where a represents 
the angle between the tangent vector t and the helical axis. This choice 



17 



ensures that the 2;-axis is parallel to the direction of the helical axis. The 
centerline of the helical rod is a space curve x(s) = {x{s),y{s), z{s)) that 
is obtained by integration of the tangent vector t(s). We have placed the 
initial position of the rod at a;(0) = 1/A;i(0), y{0) = and z{0) = (in scaled 
units), where /ci(0) is the curvature of the projection of the helical centerline 
onto the plane perpendicular to the helical axis. It is possible to show [HO] 
that 

xiO) = ^= ^^(Q^^'^Q) (50) 
^ ' h{0) Mi + A;2 (0)52(0) ■ ^ ^ 

kp{0) and M3 are free parameters that have been chosen so that the helical 
solutions displayed in Fig.|21have the same angle a. The parameters ^^^(O) = 
0.24 and M3 = 0.05 give x(0) ~ 4 for the helical solutions displayed in Figs.|2K 
and 12135 and the parameters kp{0) = 0.22 and M3 = 0.05 give a;(0) ~ 4.36 
for the helical solution displayed in the Fig. Efc. 



5 Conclusions 

The existence of circular (planar ring) and helical configurations for a rod 
with non-constant bending and twisting stiffnesses has been investigated 
within the framework of the Kirchhoff rod model. We are interested in find- 
ing analytical solutions, and since the static Kirchhoff equations consist of 
nine non-linear PDEs, we derived from them a set of differential equations 
(Eqs. (HHHini)) for the curvature, kp, and torsion, Tp, of the rod centerline, 
and then applied the geometric conditions to obtain circular and helical so- 
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lutions. 

We have shown that the planar ring {kp = Constant) is an equihbrium 
solution of the Kirchhoff equations for an inhomogeneous rod which possesses 
the bending stiffness, B{s), satisfying Eq. (j^ . and fs^s) satisfying Eq. (jHUj) . 
The planar ring is the only planar solution for twisted rods (M3 7^ 0). Non- 
twisted rods (M3 = 0) may present other types of planar solutions with kp{s) 
(non- const ant) satisfying Eq. (j2Hl)- 

We also have shown that the circular helix is a solution of the static 
Kirchhoff equations only if the rod bending coefficient B = Constant (see 
Eq. ()43|) ). An inhomogeneous rod has B ^ Constant, therefore it cannot 
have a circular helix as equilibrium solution. Nevertheless, it can exhibit gen- 
eralized helical structures that satisfy the Lancret's theorem, Eq. (jH^ . In the 
Kirchhoff model, these generalized helical solutions can be obtained solving 
Eq. (jHHjl for tf{s). Figure El displays three examples of helical structures 
that satisfy the static Kirchhoff equations. 

May an intrinsically straight homogeneous elastic rod present other types 
of generalized helical equilibrium solutions ? A homogeneous rod has C 
constant, implying = Constant (from Eq. ()13p). It has been proved that 
^' = for a helical solution of a homogeneous rod (see reference jSHI), and 
Eq. fll7|) shows that the torsion Tp = k^ = Constant. In order to satisfy 
the Lancret's theorem (Eq. (jH^ ). the curvature kp of this helical solution 
must also be constant. Since the hehcal structure whose curvature, kp, and 
torsion, Tp, are constant, is a circular helix, the only type of helical solution 
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for an intrinsically straight homogeneous rod is the circular helix. 

Some motivations for this work are related to defects jHB] and distor- 
tions j2Zl in biological molecules. Also, the stiffness of the DNA molecule 
has been proved to be sequence- dependent jS] . These defects, distortions and 
sequence-dependent elastic properties could be modeled as inhomogeneities 
along a continuous elastic rod. 

Haijun and Zhong-can j34!] have used the Goriely and Tabor stability 
analysis method [^ 1^ to explain the onset of instability in the phenomenon 
of kink formation in small closed DNA molecules reported by Han et al |38| 
1^ . In the Han et al experiment [32], the DNA rings are immersed in a 
solution containing Zn^"*" or Mg^+ ions. They found that DNA rings become 
kinked when the concentration of Zn^'^ ions is above a critical value, with 
or without the presence of Mg^"*" which alone did not produce the kinking 
phenomenon. Haijun and Zhong-can proposed that the binding of Zn^+ ions 
to DNA basepairs enhance the intrinsic curvature of the DNA destabilizing 
it. The Mg^+ ion, by binding to the phosphate backbone, did not produce 
intrinsic curvature in the DNA molecules. Nevertheless, Haijun and Zhong- 
can theory did not explain how the binding of Zn^+ ions lead to a variation 
of the curvature of the closed DNA. Also, the asymmetries and positions of 
the kinks in the DNA rings have not been explained. 

Our approach can give a new insight to the phenomenon of kink formation 
of these closed DNAs. Due to the binding of the Zn^+ ion to a basepair, the 
bending stiffness of the DNA molecule should decrease locally, leading to a 
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local increase of curvature. On the other hand, the binding of the Mg^"*" ion 
to the phosphate backbone cannot change the stiffness of the molecule since 
the phosphate backbone is a very rigid structure. Our model supports this 
proposal since the equilibrium solutions represented by Eq. show that 
the curvature is inversely proportional to the bending stiffness. Also Eq. 
(13 1|) is only valid for a non-twisted rod which is the case of the DNA rings 
reported in Ref. 

Haijun and Zhong-can have shown that closed DNA with curvature above 
a critical value never forms kinks in presence of Zn^"*" or Mg^"*" ions. It 
explains why Han et al [39j frequently observed kinks in closed DNAs of 168 
basepairs (bps) long, while they rarely observed kinks in closed DNAs of 126- 
bps long, independent of the presence of Zn^+ ions in the solution containing 
the DNA rings. However, Haijun and Zhong-can approach cannot explain 
the approximate elliptical shapes of the 126-bps circles observed by Han et 
al jHni- Our results suggest that the non-circular shape of the 126-bps closed 
DNAs is a consequence of the binding of Zn^+ ions to the basepairs that 
decreases locally the stiffness of the DNA, thus increasing the curvature. 

It has been proposed jlOl El that the kinking phenomenon is an im- 
portant mechanism for wrapping DNA around nucleosomal proteins. In 1997, 
Luger et al [-ISJ reported the X-ray crystalline structure of the nucleosomal 
protein, demonstrating that the DNA is not uniformly bent around nucleo- 
some, but presents minimum and maximum curvatures at different positions. 
It is well known that the stiffness of DNA is sequence-dependent. Our results 
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show that the variation of the bending stiffness of the DNA can give rise to 
these minimum and maximum curvatures. 
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A Appendix: The differential equations for 
the curvature and torsion 

Here, we shall derive the Eqs. (fTHH^ . Substitution of Eqs. (|T3HT7j) into 
Eqs. dHHH gives: 







f[-f2ie + rF) + f3kFCOs^ = 


0, 


(51) 






f^ + fi{C + rF)-f3kFsm^ = 


0, 


(52) 






/a - /i f^F cos^ + /2 kF sin^ = 


0, 


(53) 


(Bis) 


'kp sinO' + 


~ B{s)) kF cos a^' + rF)-f2 = 


0, 


(54) 


(Bis) 


' kp cos 0' - {C{s) 


-Bis))kFsmae + rF) + fi = 


0, 


(55) 






{C{S) {^' + TF)y = 


: . 


(56) 
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First, we extract /i and /2 from Eqs. (|55|) and (j54j) . respectively: 

/i = -iB{s) kpYcos^ + [M3 kp - B{s) kpTp] sine , (57) 
/a = (5(s)M'sine+ [M^kp - B{s) kpTp] cos ^ , (58) 

where M3 = C{s) + Tp) is the torsional moment of the rod that is constant 
by Eq. (fSUj) . Differentiating fi and /a with respect to s, substituting in Eqs. 
(1^ and (jS2I), respectively, and using Eqs. (jFrf) and (|3Hj). gives the following 
equations: 

{-iB{s)kF)"-Tp[Mskp-B{s)kpTF] + f^kF}cOS^ + 

(59) 

{[Makp - B{s) kpTp]' - tf{B{s) kp)'} sin^ = , 

{iB{s)kF)" + Tp[M^kp- B{s)kFTF]- fskp} sin^ + 

(60) 

{[Mskp - B{s) kpTp]' - tf{B{s) kp)'} cos^ = . 

Multiplying Eq. (Eq. (jHUl)) by sin^ (cos,^) and then adding the resulting 
equations, we obtain the Eq. (jl8p for the curvature and torsion: 

[M-ikp- BkpTF]' -{BkFyTF = . (61) 

Multiplying Eq. (jSHI) (Eq. (jHOI)) by — cos^ (H-sin^) and then adding the 
resulting equations, we obtain the Eq. (fT^: 

{B kp)" + kp TpiMs - Btf) - hkF = ^ ■ (62) 
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Finally, the Eq. (pUj) is obtained by substituting Eqs. (jHTj) and (fHHj) in Eq. 

(i?M' + /3 = 0. (63) 
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